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Abstract . Earth-based laser ranging to artifi- 
cial satellites and to the moon is considered as 
a technique for monitoring the Earth's polar mo- 
tion and diurnal rotation. The kinematics of Earth 
rotation as related to laser ranging is outlined. 
The current status of laser ranging as regards its 
measuring capabilities is reviewed. Artificial sa- 
tellite laser ranging has recently yielded pole 
position to better than 0.02 arcseconds with 5 
days averaging as the best result. In recent years 
single-station lunar laser ranging has produced 
UTO-values to better than 1 msec. T;i^ relative 
merits of artificial satellite and lunar laser 
ranging are pointed out. It appears that multi- 
station combined artificial satellite, and lunar 
laser ranging is likely to ultimately meet a 0.002 
arcseconds in pole position and 0,1 msec in UTl 
daily precision requirement. 

Introduction 

Although the most intriguing implications of 
the Earth orientation phenomenon lie in its dyna- 
mics and the geophysical effects involved, the 
present contribution will focus on its kinematics. 
This is done in relation to the use of Earth-bas- 
ed laser ranging to the moon and to artificial sa- 
tellites as a technique to monitor polar motion 
and diurnal rotation. It should be pointed out, 
however, that such monitoring can only be perform- 
ed efficiently if a dedicated program of measure- 
ment is based on a profound qualitative under- 
standing of the phenomenon. 

Earth orientation is involved where measure- 
ments connect positions of objects, related to an 
extra-terrestrial frame of reference , to terres- 
trial objects. The geometric relationships esta- 
blished by the measurements depend on the orienta- 
tion of a conventional Earth-fixed frame of refe- 
rence with respect to the extra-terrestrial one. 
Such situation arises in space geodesy in general 
and in satellite geodesy in particular, the moon 
conveniently being considered as a satellite. Un- 
til the advent of precise tracking means, inclu- 
ding laser ranging, the time-dependent orientation 
of the Earth was numerically modelled in terms of 
the current theories of precession and' nutation 
and the orientation parameters as provided by the 
Bureau International de I'Heure (BIH); see Veis 
[1963]. Noting the inherent precision of modern 
tracking techniques and some of their further ad- 
vantageous characteristics it was realized that 
certain constituents of Earth orientation should 
be rather modelled in terms of unknown solution 
parameters in order to exploit the qualities of 
the measurement data. One of these techniques is 
laser ranging to the moon and to artificial sa- 
tellites . 

Before reviewing the current status of this 
technique and its foreseen development the geome- 
try and kinematics of Earth orientation as per- 
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taining to laser ranging will be presented. The 
potentialities to monitor Earth orientation by la- 
ser ranging will be assessed and results obtained 
from both lunar and artificial satellite ranging 
will be quoted. These concern the diurnal rotation 
of the Earth and polar motion. 

Earth orientation as related to laser ranging 

A meaningful approach to Earth orientation re- 
quires the operational definition of at least 
two frames of reference: one to which the motion 
of the Earth is sufficiently well modelled (the 
rectangular Cartesian conventioned "Earth-fixed" 
frame: x^, Xg) and another rectangular Carte- 
sian frame z , z^) with respect to which 

the orientation is to be monitored. The defini- 
tion of reference frames in a contemporary context 
is not always straightforward as demonstrated by 
Kolaczek and Weiffenbach [ 1974] . Of particular re- 
levance is the problem of fixing a reference 
frame to a deformable Earth. The problem of Earth 
orientation is essentially the problem of finding 
the relative orientation of two reference frames, 
the components of unit vectors referred to both 
reference frames being given as time-dependent 
quantities. The basic problem is thus to find both 
sets of components , in this case from the laser 
ranging . 

In the classical optical stellar approach to 
monitor Earth orientation, the x-system is defin- 
ed by a set of conventionally adopted latitu- 
des <|> and longitudes X of participating observa- 
tories, these observatories measuring the time-de- 
pendent direction components of the observatory 
verticals with respect to a fundamental stellar 
reference frame adopted as z-frame. This z-frame 
is in fact defined by conventional positions and 
proper motions of fundamental stars and consider- 
ed inertial. Satellite ranging is, however, what 
Newton [19741 called, a "blind" technique and the 
relation to the fundamental z-frame is less 
straightforward. The main complication of the sa- 
tellite approach is that the orbits of the' sat 
tellites with respect to which the Earth's orien- 
tation is to be monitored by Earth-based measure- 
ments , have to be monitored themselves by means 
of such measurements. This holds in particular 
for artificial satellites somewhat less for the 
moon. The critical issue here is. that satellite or- 
bits referenced to the z-frame cannot be deter- 
mined from Earth-based tracking without the S 
priori involvement of the Earth's orientation. 
This involvement is however restricted so that 
Earth-based determined satellite orbits are at 
least in certain respects independent of S priori 
knowledge about the Earth's orientation, and 
Earth-based tracking data may thus contain signa- 
tures of Earth-orientation. 

Of crucial importance to any practical approach 
to Earth orientation is the operationality of re- 
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Fig. 1. Definition of reference frames by angu- 
lar Keplerian orbital elements u, i and 


ference frame definitions. The orientation of a 
right-handed rectangular Cartesian frame is uni- 
quely defined by two orthogonal unit vectors 
and Any frame obtained by rotation from this 
basin'^ frame, can then be specified by angles u , 
i* and ; see Figure 1. So is a quasi Earth-fix- 
ed x^, X -frame, in which laser tracking lo- 
cations and their modelled motions due to tides, 
tectonics etc, are specified. Likewise any y^, y_ 
y_-frame with respect to which the orientation or 
the x-frame is to be described: u^, i^, 

Given u^ , i^, , the differences 

Au = u - u-' 

2i = i^ - i^ 

Art = 


uniquely describe the rotation from the y-frame to 
the x-frame. 

Now u, i and respectively can be identified 
with the angular Keplerian orbital elements of a 
satellite: argument of latitude oj + f, inclina- 
tion and argument of the ascending node, f being 
the true anomaly and (o the argument of perigee. 
Thus the use of satellites to monitor Earth orien- 
tation seems rather obvious if u^, i^, and Au, 
Ai, An can to sufficient accuracy be determined 
as time varying quantities. The procedure to ob- 
tain such quantities is a complicated one, to be 
summarized as follows, (see Figure 2). 

Osculating' Keplerian elements of the satellite 
(t), i^(t),-o)^(t), a^(t), e^(t) agd f’^(t) are 
obtained from a state vector x(t); Jc(t), determin- 
ed in the x-frame from tracking data, e.g. laser 


ranging, in a purely kinematic way. 

An inertial y-frame can be define^ by adopting 
values u^(0), i^(0), n^(0) for u^, i , at a se- 



Fig. 2. Orientation of an Earth-fixed x-frame 
relative to an inertial y-frame by means of per- 
turbed angular Keplerian orbital elements . is 

the field of perturbing forces as described on 
the y-frame, a is the satellite orbitte semi-ma- 
jor axis, e its eccentricity and f its true 
anomaly . 


lected epoch t^. Together with a’^CO), e*(0) and 
f (0) these provide the initial state vector at t^ 
of the satellite in the y-frame. 

The two sets of angular elements 

u’^CO), i^'co), n'^(o) 

u^(0), i^(0), R^CO) 

together define the relative orientation of the 
X- and y-frames at t^: ^ = R (0) • S. 

To track the time-dependeX? orientation R (t) 
of the x-frame relative to the y-frame, anguYar 
elements 

u^(t), i^(t), fi^(t) 
are required in addition to 

U^t), i'^(t), ft^t) 

as obtained by measurement in the x-frame. Ele- 
ments 

u^(t), i^(t), R^(t) 
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can be extrapolated from the initial state vec- 
tor at tQ if the_j_complete field of forces acting 
in the y-frame (F^) is assumed. In this extrapo- 
lation the y-frame may be considered inertial. It 
should be noted however that this inertial frame 
does in general not coincide with the fundamental 
stellar z-frame of reference, which was also con- 
sidered inertial, but which deviates from the y- 
frarae by a time-invariant relative orientation 



the elements of which remain as yet unspecified. 

The force field as referred to the y-frame will 
consist of two classes of contributions: 

- forces depending on the time-dependent relative 
orientation R (t) of the x- and y-frames , e.g. 
non-central terrestrial gravitation including 
solid-Earth and ocean tides; 

- forces independent of this orientation, e.g. 
solar radiation pressure and luni -solar gravi- 
tation. 

The first class of forces poses a theoretical 
complication because the relative orientation 
R (t) to be derived has to be known in advance 
iX^order to extrapolate the initial angular ele- 
ments 

U^(0), i^(0), Q^(0) at tp 

fo obtain instantaneous values : 

u^(t), i^(t), fi^(t) at t. 

It should be noted that R -dependent inertial 
forces acting in the x-frXme do not interfere , 
provided^the procedure to obtain angular elements 
u (t), i (t), U (t) is indeed a purely kinematic 
one. 

To obtain the R -dependent force contribution 
with respect to th^^y-frame poses a comnlicationi 
leading to the introduction of what Lambeck 
[1971] called "dynamic perturbations". Although 
these perturbations may not be entirely negligible 
with future precise laser ranging to artificial 
satellites, this complication is disregarded here. 
In doing so a unique latent opportunity to verify 
the coincidence of the adopted x^-axis of maximum 
inertia [Melchior, 1972] is likewise disregarded. 
On the other hand it seems unlikely that even ad- 
vanced tracking precision would permit this veri- 
fication in a foreseeable future [Gaposchkin, 

1972; Kolaczek and Weiffenbach, 1974]. Neverthe- 
less in a detailed discussion the formal non- 
coincidence should not be overlooked: 


will be close to it ; 

c. the Earth's diurnal rotation as measured by the 
siderial angle 0 ; 

d. the small deviation termed "sway" of the Earth's 
instantaneous rotation axis from the direction 
or the total angular momentum vector; 

McClure [1973] found that for a deformable 
Earth this deviation can approach 0.01 second 
of arc. 

e. luni-solar forced nutation and precession of 
the Earth's total angular momentum vector; 

f. the relative orientation of the inertial y-frame 
to the quasi-inertial fundamental stellar z- 

f rame . 

Of these only items (b) and (c) pertain directly 
to the present subject. Ignoring for the sake of 
simplicity of the present treatment the possible 
deviations between the ^-axis and the CIO and 
between the CIO and the x^-axis and also sway, we 
can summarize: 



R, N and P being rotation matrices describing com- 
bined polar motion and diurnal rotation, forced nu- 
tation and precession respectively, in a form as 
presented by Veis [1963]. Writing 



P ' may be considered as describing precession with 
respect to the arbitrarily defined inertial y- 
frame, leading to the conclusion that satellite 
techniques yield the compound rotation 

y = P' . N ,R . 5 , 

precession P' relative to the y-frame. To comply 
in practice with convention and to conveniently 
separate R from precession and nutation, when 
assuming the latter two, the initial angular or- 
bital elements 

u^(0), i^(0), n^(o) 

are selected so, that the y-frame coincides with 
the z-frame , although because of incomplete know- 
ledge about R(0) at epoch, this can be realized 
only approximately, even if N(0) and P(0) are 
assumed. If realized, R is unity and P' = P, so 
that: y" 

P . N , R = R 

zx 

and consequently: 


x5 


the C-system being the axis-of-figure (maximum 
inertia) system. Hence: 




■ 5 


R ,, thus defined includes: 

Zi 

a. 'the unknown small and virtually invariable de- 

viation of the adopted quasi Earth-fixed x- 
frame from the axis-of-figure g-frame; 

b. polar motion relative to the x^^-axis, which 
does not necessarily coincide with the CIO, but 


R(t) = N* 


R 


If not fully realized, R(t), thus derived, will be 
biased by an unknown, but contstant relative 
orientation. 

A sequence of R(t) for a sequence of instants 
t will provide the time-dependent orientation of 
the Earth's rotation axis with respect to the 

Earth-fixed x-frame (:ji , X ) and the siderial 
~D O 

angle 0, see Figures 3*^ and 4. Finally, to comply 
with IPMS and BIH convention, pole position is 
given as x and y . 

D ^ "D XXX 

In practice [Lambeck, 1971] u , i , may be 
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Fig. 3. Geometry of Earth orientation with 
respect to stellar z-frame, the orientation to 
be obtained from satellite orbits. X and <p 
specify the direction of the Earth ' s^rotatiBn 
axis relative to an Earth-fixed x-frarae. 


transformed to u^', as referred to an 

intermediate quasi inertial z' -frame, which appro- 
ximates the z-frame: 


t = D(t) . 


This is performed by approximating R(t) by R(t), 
replacing 0 by an approximate value 6 and equating 


X and y to zero: 
P P 


cos 9 -sin 6 0 


R(t) 


sin 9 cos 9 0 

0 0 1 


Solving then for D(t), rather than for R(t) direct- 
ly, yields'!' , X and A9 = 0-9; see Figure 5. 

This geneBal ^ approach is valid for both ar- 
tificial satellites and the moon although it is 
primarily tuned in to the use of artificial sa- 
tellites, not restricted however to laser ranging. 
Before pointing out special features of the lunar 
laser ranging case some general conclusions may be 
drawn; 

.- because the x-frame does not necessarily coin- 
cide with the conventional terrestrial reference 
frame as' defined by the CIO and the conventional 
zero meridian of Greenwich, there may appear 
constant biases in all three determined rotatio- 



Fig. 4. From the relative orientation of Earth- 
fixed X- and stellar z-frames to pole position 
(x ,y ) and siderial angle 0. 


loP.N.R. X 


I'cxP.N.R.x 


R= D.R = D. 


co*0 ~sinO 0 
s i n 9 coo 9 0 

0 0 1 



Fig. 5. Pole position (x ,y ) and siderial angle 
(0) from the relative oriintStion of the quasi- 
inertial z' -frame relative to the inertial z- 
frame. 0 is an approximated siderial angle. 
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nal parameters Xp, yp, 9 as derived in the 
classical optical stellar way; 

- due to the well-known longitude-ambiguity in ar- 
tificial satellite orbits such bias is so likely 
in 0 that artificial satellite techniques are 
expected not to provide UTl, but rather the ro- 
tation rate lo of the Earth; 

- additional not necessarily constant biases may 
be caused by erroneously modelled lunar-solar 
precession and mutation (P-and H-matrices); 

- a detailed force model and a scrupulous applica- 
tion of that model are required. 

The moon is just another satellite, kinematically 
differing from the others by: 

- its larger mean orbital radius; 

- its relatedly longer period of revolution, in 
fact much longer than that of the Earth's diur- 
anl rotation; 

- its more stable center of mass orbit , less effect- 
ed by the Earth's graviational potential and 
other terrestrial perturbing forces; 

- its substantial size, which requires the seleno- 
d.etic coordinates of the retroreflector sites 
and the lunar rotation to be modelled or solved 
for; 

- there is only one moon, hence one lunar center- 
of-mass orbit but the moon carries several geo- 
metrically distinct retroref lectors. 

These differences entail marked differences con- 
cerning tracking-operational and data analysis 
aspects. These are indicated in the next two para- 
graphs respectively. 

The status of laser ranging 

The technique itself is assumed known and will 
be reviewed only as regards its main character- 
istics pertaining to Earth orientation determina- 
tion. Detailed and up to date technical informa- 
tion on availabe instrumentation and current deve- 
lopments can be obtained from Pearlman and Hamal 
[1978], 

Relevant system characteristics are: 

- maximum range capacity; 

- ranging precision; 

- day-light ranging capability; 

- repetition rate. 

Speaking in general terms these characteristics 
are not mutually independent. 

In the present context, systems can as regards 
their maximum range capibility be classified in 
broadly three catagories : 

- those of lunar ranging capability; 

- those able to range at Lageos; 

- those restricted to closer artificial satellites, 
like Beacon Explorer-C, Starlette and Geos-3, 

It should be noted that the operational aspects 
of lunar laser ranging (LLR) are quite different 
from those of artificial satellite laser ranging 
(SLR), and a system of lunar range capability will 
not necessarily be able to range at artificial sa- 
tellites. 

The most advanced SLR systems attain between 5 
and 10 cm single shot precision, expressed as a 
distance standard deviation; 2 to 3 cm is foreseen 
for the next few years. LLR "normal observation 
point" precisions (see next paragraph) reach or 
are expected to reach the same level. 


Dayrlight ranging capability is of crucial im- 
portance in order to provide a continuous record 
of short period Earth orientation phenomena, such 
resolution requiring averaging times of a fraction 
of a day. It is to be understood that such short 
averaging time is a prerequisite of future precise 
Earth rotation monitoring systems. Many of the 
operational and most of the planned SLR devices 
have day-light capability on the closer satellites , 
few on Lageos. LLR features guiding difficulties 
when ranging is attempted close to new moon. Be- 
cause of that, the only routinely operational LLR 
system (at the McDonald Observatory, Fort Davis, 
Texas) cannot effectively range within 3 days off 
new moon. 

The repetition rates of most SLR systems exceed 
0.1 pps. The McDonald LLR performs 1 pulse per 3 
seconds, but the McDonald team has made it practice 
to compress the data from each 5 to 20 minutes run 
into a single "normal point", three runs being 
attempted per day, leading to an equal number of 
daily "normal points", except around new moon. 

Of paramount importance for Earth orientation 
work, will be station siting, considering average 
atmospheric conditions and the need of global de- 
pl 03 nnent, the first requirement in view of data 
continuity, the second for "Earth-fixed" (x-) re- 
ference frame definition. 

Although only part of the available facilities 
has been used in dedicated programs of polar mot ion 
and diurnal rotation studies, encouraging prelimi- 
nary results have been obtained from both SLR and 
LLR. 

Review of work accomplished 

In the past eight years or so several types of con- 
tribution have been made to the determination of 
the Earth's diurnal rotation and polar motion: 

- theoretical modelling; 

-feasibility analysis of simulated data; 

- dedicated data taking; 

- analysis of actual data. 

Only the latter of these will be reviewed, separa- 
tely for SLR and LLR, arbitrarily, in this order. 

In the present context this review can only be 
sketchy and cannot satisfactory reflect the amount 
of effort spent by contributing individual or 
groups of investigators. 

Early attempts at NASA's Goddard Space Flight 
Center (GSFC) to detect polar motion from artifi- 
cial laser ranging were initiated in 1970 and de- 
monstrated the capability. [Smith et al, 1972a]. 
First preliminary results [Smith etal, 1972b] in- 
dicated that using the data of a single 30 cm pre- 
cision ranging station the variation in latitude 
of that station could be derived to 0.03 seconds 
of arc with a time resolution of 6 hours. Consi- 
dering the latitude of that station at Greenbelt, 

MD (39 N) and the inclination of the single Beacon 
Explorer-C satellite used (41 ) this satellite 
could be traclfed near apex and i' (briefly written 
instead of i ) be accurately determined to about 
0.001 second of arc, the accuracy of Ai=i'-i li- 
mited by that of i (briefly, instead of i ) , be- 
cause of gravitational field uncertainties. What 
could be derived from a single station was the va- 
riation of latitude, thus polar motion projected 
onto the Greenbelt-meridian; not the corresponding 
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component of polar deviation from the CIO, simply 
not because the Greenbelt-latitude is not knovm 
in the conventional BIH-frame. The method employed 
to obtain such single station polar motion results 
has been described to some detail in [Dunn et al, 
1974] and [Kolenkiewicz et al, 1977], It is a spe- 
cial version of the general approach outlined be- 
fore, the angular element analysed being the in- 
clination i, the pertinent kinematic equation rea- 
ding [Lambeck, 1971]: 

i> - i = _i() sin(X„ - X ) , 

P N p ’ 

where Xj^ stands for the longitude of the satellitefe 
ascending node. Using up to four consecutive sta- 
tion passes of the satellite within a time span 
of about 6 hours, short arc osculating inclinations 
i' were obtained (see Figure 6) and compared with** 
a reference orbit , extended over the entire period 
of the experiment. Later this so-called "max-lat" 
approach was abandoned and replaced by a more fle- 
xible approach yielding both variation of latitude 
and length of day information [Dunn et al, 1977], 
During a 3-week period in 1970 a second ranging 
station operated from Seneca, N,Y., 400 km North 
of Greenbelt, yielding also a small number of 4- 
pass arcs. This enabled two independent determina- 
tions of both solution parameters. 

The single-station approach as just described, 
has several limitations; 

- of polar motion only the meridional component 
can be measured; 

- it is difficult to separate this component from 
dynamical perturbations of i, these to be mo- 
delled over the full time span of the investiga- 
tion; 

- it is impossible to separate polar motion from 
precession and nutation; 

- the use of several satellites or interrupted 
arcs of the same satellite may cause discon- 
tinuities and inconsistencies in the pole path 
as measured. 

These can completely or to a substantial degree be 
overcome by deploying a multi-station network. 

Such network enables to determine a complete set 
of osculating angular orbital elements 
from.'.a day or less of laser tracking. A minimally 
perturbed reference orbit u (t), i^(t), fi (t) va- 
lid only for the same period of a day or less, say, 
will suffice to yield a sufficiently detailed 
R (t)- or R -record to derive the direction of 
tgg orientation axis for the observation period. 
Because of the shortness of this period incomplete 
modelling of precession (E) and nutation (N) and 
of acting forces will practically not interfere. 

To be more precise as concerns precession and nu- 
tation an error in P(0) and N(0) at t would cause 
a biased reference orbit, but because°of the short- 
ness of the tracking period this bias will remain 
constant throughout and persist as such in P*(t) 
and N*(t), 

Polar motion results from multi-station laser 
tracking of artificial satellites were reported 
recently from two sources by Smith et al [1978b] 
and: Schutz et al ri97Sb]-. These results announced 
a break-through in polar motion determination from 
laser ranging. 

Smith et al [1978b] reported preliminary results 
obtained from Lageos tracking at a total of seven 


GSFC- and Smithsonian Astrophysical Observatory 
(SAO)-stations, four of which in the U.S.A., two 
in South-America and one in Australia. The four 
GSFC-stations in the U.S.A. claim a 10 cm single 
shot ranging precision, the SAO-stations in South- 
America and Australia 1 m. First a consistent set 
of station coordinates was obtained from this data 
to define an x-frame of reference [Smith et al, 
1978a]. Subsequently an iterative procedure of 
fitting 5-day arcs to three interlinking 30-day 
reference orbits yielded average pole positions, 
relative to the adopted x-frame, for each of the 
eighteen 5-day periods of October, November and 
December 1976. The formal standard error in the 
X -component of polar position is 0.003 arcseconds, 
0?002 arcseconds in y . The authors believe however 
that a precision betwien 0.01 and 0.02 arcseconds 
is more realistic. From this experience and pre- 
vious simulations [Kolenkiewics et al, 1977] it 
was concluded that with more stations ranging La- 
geos at the 10 cm level and improved modelling, in 
particular of solid-Earth and ocean tidal effects, 
ultimately daily 5-cm pole position and UTl values 
to about 0.2 msec consistent over 2 to 3 months 
are feasible. Recently GSFC issued a first bulletin 
of preliminary 5-day pole positions as obtained 
from Lageos ranging; the period covered is May-De- 
cember 1976. Standard deviations given there range 
on the average from about 0.01 to about 0.03 arc- 
seconds . 

Schutz et al [1978b] reported pole positions ob- 
tained from laser tracking of Geos-3, spanning a 
roughly one-month period in early 1976. Data was 
used from three GSFC 5 to 10 cm single shot pre- 
cision systems located at Greenbelt, MD, Bermuda 
and Grand Turk. The analysis was based on a pre- 
viously obtained set of station coordinates de- 



Fig. 6. A single station P (rotated to respec- 
tive positions P^, P^, Pg and P ) is able to 
track four consecutive passes of an artificial 
satellite near apex and thus derive its orbital 
inclination. 
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fining an x-system and the GEM-10 geopotential mo- 
del. The data span was devided into sixteen 2-day 
arcs . Special care was taken to model upper at- 
mosphere effects. Weighted least squares straight 
line fits to the determined x and y values yield- 
ed estimated x and y standard deviStions smaller 
than 0.07 arcsiconds .^The results indicate that a 
determination of both pole position components is 
possible from a regional station network, that as 
expected, the y -component is more accurately re- 
covered than thE x -component and that additional 
model improvementsPmay be required for further 
analysis of Geos-3 data. The method used has been 
documented in more detail in connection with po- 
lar motion analysis of earlier Geos-3 data [Schutz 
et al, 1978a}. Although at first glance not im- 
pressive the authors consider the results from the 
1976 data analysis of particular significance be- 
cause this data was not used in the GEM- 10 geopo- 
tential and station coordinates solution. 

The LLR experiment has been outlined by Bender 
et al [1973]. A total of five retroflector arrays 
has been deposited on widely dispersed sites on 
the moon by the US and USSR space programs and 
subsequently laser ranged. Although returns have 
been announced by nine stations in five countries 
[Silverberg, 1978], LLR on a routine basis has in 
essence so far been a single-station operation from 
the McDonald Observatory. This restricted the 
yield of practical results as regards Earth orien- 
tation, like this was restricted in single-station 
SLR. Nevertheless, like single-station SLR demon- 
strated a capability to measure polar motion, so 
demonstrated the single-station LLR a capability 
to measure UTO as pertaining to that station . 

The lunar orbit differs most markedly from ar- 
tificial satellite orbits in its much larger semi- 
major axis, which entails a 200 to 400 times longer 
period of revolution and a superior orbital co- 
varage from a single observing station. These fea- 
tures have consequencies for both the data acqui- 
sition and the data analysis . Apart from a gap 
around new moon as caused by guiding problems which 
seem difficult to overcome , a geographically well 
selected LLR station can observe the moon daily, 
more or less evenly distributed over a full revo- 
lution. As compared to SLR the longer period and 
the high degree of predictability of the lunar or- 
bit entail distinct differences in observation da- 
ta handling. A LLR "normal point" covers about 5- 
15 minutes of ranging, a portion of a lunation too 
short to construct orbital elements. Such normal 
point may be regarded as an equivalent photon re- 
turn representing the observation run. Over the 
past few years an average of 25-30 normal points 
have been produced at McDonald per lunation [Mul- 
holland, 1978]. When solving for Earth orientation 
a normal point is fitted to a retroreflector re- 
ference orbit , obtained by taking into account the 
selenodetic retroreflector coordinates and the 
physical librations of the moon. The range resi- 
duals are then analysed in an attempt to recover 
Earth orientation signatures. This procedure 
follows the general pattern as outlined before 
and seems to deviate from SLR data analyses to 
recover Earth orientation, only in so far dictat- 
ed by the kinematics of the lunar orbit and the 
ensuing data acquisition strategy. 

Range residuals have been analysed by several 


authors to study various aspects of the Earth-moon 
system, including the Earth's diurnal rotation in 
terms of UTO; Stolz et al [1976], Harris and 
Williams [1977], Shelus et al [1977], King et al 
[1978], Calame [1978]. Single-station operation 
implies that as regards the direction of the ro- 
tation axis, only the projection of its change 
(polar motion) with respect to an Earth-fixed x- 
frame, onto the station's meridian can be measur- 
ed. This, in turn, implies that not UTl, but only 
UTO can be derived, since the transverse component 
of pole position, required to correct UTO in order 
to obtain UTl, cannot be determined from a single 
station. Therefore Earth orientation results from 
LLR have up till now only been obtained in terms 
of UTO-corrections valid for McDonald. This never- 
theless is a unique contribution, because it is a 
field in which reliable results from SLR are not 
expected in the near future , and if so , these will 
measure the rotation rate lo, rather than the angu- 
lar position of the Earth-fixed frame in terms of 
the siderial angle as obtained from UTl. 

First UTO-values derived from LLR were reported 
by Stolz et al [1976]. Using diurnal variation of 
range residuals they obtained, scattered over a 
roughly 5 years' time span, 194 single day values 
with a medium standard deviation of 0.7 msec if 
allowance is made for the uncertainty of the 
current lunar ephemeris . Only days with data well 
distributed in time were included in the analysis. 

Rather than selecting the best isolated obser- 
vation days of a data span, Shelus et al [1977] 
obtained UTO results for a complete lunation on a 
daily basis, accepting all observations, as one 
would have to do in a routinely operating Earth 
rotation monitoring service . The results were ob- 
viously less precise than those of Stolz et al 
[1976], the standard deviation relative to the 
corresponding BIH data amounting to 2.6 msec in 
the most realistic case of analysis, but they may 
provide a more realistic measure of what can be ex- 
pected on a continuing near real-time basis from a 
single station operating under normal conditions. 
The results obtained by Stolz et al [1976], on the 
other hand, demonstrate what can be done on a dai- 
ly averaging basis under the best conditions. It 
should possibly not be overlooked that although 
Shelus et al [1977] did not select the best days 
of data within a lunation, they seem to have se- 
lected the best LLR lunation at McDonald . 

King et al [1978] analysed the normal points on 
four retroreflectors taken at McDonald between Oc- 
tober 1970 and November 1975 and solved for correc- 
tions to UTO in terms of 126 "tabular points" de- 
fining a continuous, piecewise linear function of 
time spanning the 5-year interval. The solution 
was a phased one for a total of 166 parameters , 
only the UTO parameters being obtained in the se- 
cond phase . After removing a constant difference , 
the standard deviation of the tabular points as 
compared to BIH values is 2.1 msec. The standard 
deviation of a, comparison of the tabular points 
with the results obtained by Stolz et al [1976], 
who analysed most of . the same data, is 0.8 msec. 

Calame [1978] proceeded along several lines, 
varying the selection of parameters solved for si- 
multaneously with UTO, data selection criteria and 
averaging time (1, 2 or 5 days). Using data 
through January 1978, the estimated UTO standard 
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deviation ranges from about 1 to 3 msec. 

A major break-through as regards Earth orienta- 
tion from LLR is expected from the deployment of 
the proposed multi-station EROLD-network [e.g. 
Mulholland and Calame, 19781. Such network would 
enable the separation of all three components of 
Earth orientation at a level of precision assessed 
by Stolz and harden [19.77]: better than measuring 
accuracy with an averaging time of two days. The 
BIH's involvement in EROLD has been outlined by 
Calame et al [1976]. 

Summarizing one could say that both SLR and LLR 
have already demonstrated some of their capabili- 
ties and defaults in single-station operation to 
measure Earth orientation. Multi-station results 
are available from SLR. Multi-station LLR is ex- 
pected in the near future. Up till now SLR and LLR 
have been complementary in that SLR provided main- 
ly polar motion, while LLR yielded information on 
diurnal rotation. In multi-station operation one 
would expect both SLR and LLR to provide complete 
three-parameter Earth-orientation, comprising both 
polar motion and diurnal rotation. Short averaging 
times are more likely to be achieved with SLR. On 
the other hand, the lunar orbit offers a more 
stable long term reference, in particular for the 
determination of UTl and ui. Considering this, SLR 
and LLR should be considered complementary tech- 
niques, rather than competitive. It is important 
to note here that at least some of the planned 
LLR stations will also have Lageos-capability. Al- 
ternatively Silverberg [1978] envisions SLR sta- 
tions with occasional LLR capability. In such in- 


tegrated arrangements the SLR could be used to 
track polar motion on a daily basis, the LLR from 
one or more of the stations being then used to 
measure UTl. applying the SLR polar coordinates to 
correct UTO as measured. Ultimately Earth orien- 
tation to equivalent measurement precision or 
better with averaging times of one day or shorter 
is expected from integrated LR. 

Concluding remarks 

Figure 7 depicts types of Earth orientation 
results obtained or expected from laser ranging, 
classified according to modes of operation. As 
pointed out in the preceding paragraphs, only 
multi-station operation will be able to provide 
both polar motion and UTl results. Ultimately 
both SLR and LLR should be able to do so indepen- 
dently. Up till, the present time, only SLR has 
demonstrated multi-station operation and yielded 
polar motion in this mode. Once in multi-station 
operation, as foreseen in the EROLD campaign, LLR 
is more likely than SLR to provide long-term con- 
sistency of results in particular as regards UTl. 
On the other hand SLR may offer shorter averaging 
times. Considering moreover practical operational 
constraints , there is strong tendency to foresee 
SLR and LLR as operating in a complementary ra- 
ther than in a competitive way when monitoring 
Earth orientation, allowing both techniques to 
contribute on their strong points. Instruments 
having both SLR and LLR capability seem advanta- 
geous in this respect. A future Earth orientation 
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Fig. 7. Classification of Earth orientation results. Solid lines indicate demonstrated capabilities, 
dashed lines indicate potentialities. 
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service is supposed to provide pole position to 
0.002 arcsec and UTl to 0.1 msec with averaging 
times of one day or shorter. When extrapolating 
current experience into the future by means of 
simulation [e.g.. Stolz and harden, 1977; Smith et 
al, 1978b] it seems that the above requirements 
can be ultimately met by multi-station combined 
SLR and LLR. Additional requirements to be con- 
sidered for a service bureau type of operation 
are continuity of results, quasi real time 
availability of results and cost-effectiveness. 

Acknowledgement . The author is indebted to Dr. 
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